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Abstract 



We discuss how a possible violation of the combined symmetry CPT in the 
B meson system can be investigated at the LHC We show how a tagged and 
an untagged analysis of the decay modes of both B d and B s mesons can lead 
not only to a possible detection of a CPT-violating new physics but also to an 
understanding of its precise nature. The implication of CPT violation to a large 
mixing phase in the B s system is also discussed. 
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1 Introduction 



The combined symmetry CPT is supposed to be an exact symmetry of any local axiomatic 
quantum field theory. This is indeed supported by the experiments: all possible tests so far 
PQ have yielded negative results, consistent with no CPT violation. Why then should we 
be interested in the possibility of CPT violation in the B system? There are three main 
reasons: first, any symmetry which is supposed to be exact ought to be questioned and 
investigated, and we may get a surprise, just like the discovery of CP violation; second, it 
is not obvious that CPT will still be an exact symmetry in the bound state of quarks and 
antiquarks, where the asymptotic states are not uniquely defined j2]; third, there may be 
some nonlocal and nonrenormalisable string-theoretic effects at the Planck scale which have 
a ramification at the weak scale through the effective Hamiltonian [3]. Moreover, this effect 
can very well be flavour-sensitive, and so the constraints obtained from the K system [4] may 
not be applicable to the B systems. A comprehensive study of CPT violation in the neutral 
K meson system, with a formulation that is closely analogous to that in the B system, may 
be found in [5]. 

There are already some investigations on CPT violation in B systems. Datta et al. [6] 
have shown how CPT violation can lead to a significant lifetime difference AT/T in the 
generic M°-M° system, where M° = K°,B°, or B s . It was discussed in [7] how direct 
CP asymmetries and semileptonic decays can act as a probe of CPT violation. Signatures 
of CPT violation in non-CP eigenstate channels was discussed in [8]. The role of dilepton 
asymmetry as a test of CPT violation and possible discrimination from AB = —AQ processes 
were investigated in j9] . The BaBar experiment at SLAC has tried to look for CPT violation 
in the diurnal variations of CP- violating observables and set some limits [TO] . 

Right now, there is no signature of CPT violation, or for that matter any type of new 
physics, in the width difference of B° — B° and decay channels of Bd 0- The width difference 
for the Bd system, Ar^, is too small yet to be measured experimentally, and the bound is 
compatible with the Standard Model (SM). On the other hand, it is expected that the width 
difference Ar s would be significant for the B s system, but at the same time we know that 
the theoretical uncertainties are quite significant [TT]. If there is some new physics (NP) that 
does not contribute to the absorptive part of the B° s —B s ° box, the width difference can only 
go down [12], while there are models where this conclusion may not be true [13]. To add to 
this murky situation, the CP- violating phase /3 S , which is expected to be very small from the 
CKM paradigm, has been measured [H] to be large, compatible with the SM expectations 
only at the 2.1a level. The situation is interesting: there is hint of some NP, but we are yet 
to be certain of its exact nature, not to mention whether it exists at all. 

1 We use B° and B Q to indicate the flavour eigenstates, Bd as a generic symbol for both of them, and 
similarly for B s . The symbol B q will mean either a Bd or a B s . 



In this situation, let us try to see what we can expect at the LHC, where the B s meson, 
along with the B d , will be copiously produced. We are helped by the fact that the time 
resolution in ATLAS and CMS are of the order of 40 fs, so one can track the time evolution of 
even the rapidly oscillating B s . Thus, we expect excellent tagged and untagged measurements 
of both B d and B s mesons. It is best to focus upon the single-amplitude observables: B d — > 
J/ipKs and B s — > J/ip<f) or B s — > D+D~ For the J/ip<f) mode, one has to perform the 
angular analysis and untangle the CP-even and CP-odd channels. 

In this paper, we will discuss how one can detect the presence of a CPT violating new 
physics from the tagged and untagged measurements of the decay. We will confine our 
discussion to the case where CPT violation is small compared to the SM amplitude, just to 
make the results more transparent. The conclusions do not change qualitatively if the CPT 
violation is large, which, we must say, is a far-off possibility based on the data from the 
other experiments [10] . We will also show how the nature of the CPT violating term can be 
probed through these measurements. 

In Section 2, we mention the relevant expressions, and introduce CPT violation, with 
relevant expressions, in Section 3. The analysis for both Bd and B s systems is performed in 
Section 4, while we summarise and conclude in Section 5. 



2 Basic Formalism 



Let us introduce CPT violation in the Hamiltonian matrix through the parameter 5 which 
can potentially be complex: 
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where 5' is defined by 



lb' 



VH12H21 

Solving the eigenvalue equation of Ai, we get, 
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2 They are not strictly single- channel as there is a penguin process whose dominant part has the same 
phase as the leading Cabibbo-allowed tree process, but on the other hand these channels are easy to measure, 
and the penguin pollution is quite small and well under control. 
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Hence, corresponding eigenvectors or the mass eigenstates are defined as 

\B H ) = Pl \B°) + qi \B°) , 



The normalisation satisfies 
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Let us define, 
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The convention of [TU] leads to z = 5/2, where Zq is a measure of CPT violation as used 
in [10]. The limits imply that \zq\ 1. Even if the origin of CPT violation is something 
different, it is not unrealistic to assume \5\ 1. 

One could even relax the assumption of H 2 \ = H* 2 - However, there are two points that 
one must note. First, the effect of expressing H\ 2 = h\ 2 + 8, H 2 i = h\ 2 — 8 appears as 8 2 
in yJHi 2 H 2 i, the relevant expression in eq. (1), and can be neglected if we assume 8 to be 
small. The second point, which is more important, is that CPT conservation constrains only 
the diagonal elements and puts no constraint whatsoever on the off-diagonal elements. It 
has been shown in [5J [7] that H\ 2 ^ leads to T violation, and only Hn ^ H 22 leads to 
unambiguous CPT violation. Thus, we will focus on the parametrization used in eqs. (1) 
and (2) to discuss the effects of CPT violation. 

The time-dependent flavour eigenstates are given by 
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So, the decay rate of the meson B q at time t to a CP eigenstate / is given by 
r(B q (t)^f CP ) = [\f4t)\ 2 + \^fAV4t)\ 2 + 2Re^ f J4t)fm}\ A f\ 2 

r(B q (t) -+ f CP ) = [|/-WI 2 + l0 2 l 2 l/ + Wr + 2Re(OJ + W/-W)] 

where 

A, = (f\H\B q ) , A, = {f\H\B q ) . 

Also, 
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In the SM, both are equal and £/ x = £/ 2 = £/. For single-channel processes, = 1. 
Now using eq. (T7|) and eq. ([2]) one gets 
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Here, Am and Ar are defined through; 



with 
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3 Introducing CPT Violation 



If we consider a time-independent CPT violation so that 5 is a constant, there are only two 
unknowns in the picture: Re(5) and lm(5), over those in the SM. We will try to see how one 
can extract informations about them. For our analysis, let us take 5 to be complex; it will 
be straightforward to go to the simpler limiting cases where 5 is purely real or imaginary. 
For example, if the width difference Ar is much smaller than Am, the model of [10] makes 
8 completely real. 

When B q and B q are produced in equal numbers, the untagged decay rate can be defined 



The above equation is useful to fix the overall normalization. 

We assume, 5^1 and expand any function f(5) using Taylor series expansion and drop 
all the terms 0(5 n ) for n > 2. From eq. ([TBI) , eq. (fTUj) and eq. (TTHT) we will get the untagged 
decay rate for the decay B q — > /, 



as 
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using the above expression one could define the branching fraction as 
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Thus, for the B s system, where the hyperbolic functions are not negligible, we get (keep- 



ing up to first order of terms in Ar s ): 

"OO 

2 



Br[f} = ll dtV [/./] 



o 



\A 



f\ 



^{(l + le/| 2 )(l + 2^)-Im(5)Im(0)' 
+ { 2Re (0) " 5(1 " lOl 2 )Re(5) - iRe(0)((Re(5)) 2 - (Im(5)) 2 )} 



I /\ 777 

:ImW{(l-|e/| 2 ) + Re(5)Re(0)} 



(19) 



Theoretically, one can obtain the coefficients of the trigonometric and the hyperbolic 
terms by fitting the untagged decay rate. In actual cases this is a difficult task. However, 
there is one other observable which may help us. Before we go to that, let us note that the 
above expression is further simplified in the following four cases. 

• For the system: We can neglect AI^ so that the cosh term is unity and the sinh term 
is zero. Thus, there are only two time-dependent terms, cos(Amt) and sin(Amt), and 
the fitting is easier. Note that AI^ is measured to be small, so we need not consider 
the case where it is enhanced to a significant value because of the CPT violation. In 
fact, if 5 is small, AI^ is bound to be that coming from the SM, as the correction is 
proportional only to 5 2 and higher. 

• For one-amplitude processes: We can put = 1, and only one of Re(£/) and Im(£/) 
remains a free parameter 

• For 5 being either purely real or purely imaginary: The expressions are straightforward. 
For example, if 5 is purely real, there is no trigonometric dependence on the untagged 
rate. 

• Finally, for \5\ 1: We can neglect terms higher than linear in either Re(5) or Im(<5) in 
eq. f|T9|) . This is expected to be the case according to [10]. For example, the expression 
for the branching fraction for a one-amplitude process simplifies to 
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3 £/ is a SM quantity, so it is theoretically calculable, but it may also contain other new physics which is 
CPT conserving, so it is better to obtain both real and imaginary parts of £/ and check whether = 1. 



One can also tag the B mesons and define a tagged decay rate asymmetry 
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Note that (i) for Re(5)=Im(5)=0, this reverts back to the SM expression, as it should, 
and (ii) If \8\ <C 1 and AT /T <C 1 as in the B^ system, the tagged rate can measure both 
Re(<5) and Im(<5). 

For one-amplitude processes with \5\ <C 1, one may write a simplified expression: 
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It is clear from eq. ( 122]) how one can extract Re(<5) and Im(<5) by comparing the untagged 
and tagged analyses. Suppose weconsider the B s system where Ar s is non-negligible. The 
coefficient of the sinh term in IV gives Re(<5). However, there is an overall normalisation 
uncertainty given by \A f \ 2 . To remove this, one can consider a combined study of the coeffi- 
cients of sinh (— j^) and cos (Am s t) from the untagged and tagged decay rates respectively; 
their ratio allows for a clean extraction of Re(<5). On the other hand, the ratio of the co- 
efficients of cos(Am<jt) in Tjj and sin(Am s t) in IV gives a clean measurement of lm(S), as 
Im(£j) is known from the SM dynamics. A further check about the one-amplitude nature is 
provided from |Re(£/)| 2 + |Im(£/)| 2 = 1. In fact, as long as 5 is small, one can extract both 
Re(<5) and lm(5) even if ^ 1, from the coefficients of the sine, cosine, and hyperbolic sine 
terms in and IV. 



One may also define the time- dependent CPT asymmetry as 
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and the time-independent CPT asymmetry as 
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This goes to the usual CP asymmetry Acp if 5 — 0; thus, any deviation from the expected 
CP asymmetry calculated from the SM would signal new physics, but one must check all 
the boxes to pinpoint the nature of the new physics. For example, there would not be any 
change in the semileptonic CP asymmetry if the new physics is only CPT violating in nature. 

4 Analysis 

There are five a priori unknowns: Re(<5), Im(<5), Re(£/), Im(£j), and \Af\ 2 . For a one- 
amplitude process |£/| 2 = 1 and the number of unknowns reduce to four. The tagged and 
untagged decay rates, the branching fraction, and the time-independent CPT asymmetry 
would provide informations on all of these unknowns. Assuming the CPT-conserving physics 
to be purely that of the SM, one may calculate £j following the CKM picture. In the analysis 
that follows, we take £/ to be known from the SM. We would like to point out the following 
features. 

• The overall amplitude \Af\ 2 cancels in the CPT asymmetry. This, therefore, is going 
to be the observable one needs to measure most precisely. 

• It is enough to measure the coefficients of the trigonometric terms only. For the Bd 
system, Ar d is small anyway, and for the B s system, Ar s has a large theoretical 
uncertainty. 

• The analysis holds even if the process under consideration is not a one-amplitude 
process. In fact, one may check whether there is a second CPT conserving new physics 
amplitude just by looking at the extracted values of Re(£/) and Im(£y). 

• The coefficient of sin(Am g t) in the expression for the tagged decay rate T T gives 
the mixing phase in the box diagram. Thus, lm(S) may be constrained by the CP 
asymmetry measurements in the B d system. On the other hand, even those constrained 
values generate a large mixing phase for the B s system compatible with the CDF data. 



4.1 The B s system 



For the B s system, we take 

Am s = 17.77 ± 0.12PS- 1 , Ar s = 0.096 ± 0.039PS" 1 , — - = 0.147 ± 0.060 , 

J- s 

— = 1.530 ± 0.009ps , Re(^) = 0.99 , Im(£/) = -0.04 . (25) 

J- s 

In figure [H we show the variation of Aqpt with Re(<5). For our analysis, we take both 
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Figure 1: Variation of Aqpt with Re(<5) for the B s system. The three lines, from top to 
bottom, are for Im(#) = —0.1,0 and 0.1 respectively. 

|Re(5)|, |Im(5)| < 0.1, which is consistent with [TU]. The variation of Aqpt with Am s and 
Ar s is negligible, of the order of 0.2%, so we fix them to their respective central values. 
Effects of 5 in both Am s and Ar s are quadratic in 5, and hence we can use the SM values 
for them. In fact, Aqpt does not depend significantly on the choice of Im(5) either; the 
variation is less than 1%. This is due to the fact that here, |Im(£y)| |Re(£y)| and hence 
the coefficient of Re(5) is much greater than the coefficient of Im(<5) in the expression of 
Aqpt- This feature does not hold for the B d system. Note that Aqpt clearly gives the sign 
of Re (5). The small nonzero value of Aqpt for 5 = indicates the small mixing phase in the 
B® — B s ° box diagram. However, the apparent phase, i.e., the coefficient of sm(Am s t) , can 
increase with Im(<5), as can be seen from figure [21 

4.2 The B d system 

The inputs that we use for the B^ system are 

Am d = 0.507PS" 1 , AT d = , Re(^) = 0.72 , Im(£/) = 0.695 . (26) 
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Figure 2: Variation of sin(2/3 s ) with lm(5). 

This follows from the CKM expectation of sin(2/3d) = 0.695 ± 0.020. The constraint on 5 
comes from the measurement of sin(2/3^) in the b — > ccs channel: 0.668 ±0.028 pj)] 0. Again, 
we can fix Am^ at its central value. This time, due to the comparable values of Re(£/) and 
Im(£y), Acpt is sensitive to both Re(5) and Im(5). The variations are shown in figure [3] for 
three values of Im(5) and figure @] for three values of Re(<5). It turns out that A CPT is always 
positive for Re(<5), Im(<5) < 1; this is a consistency check for the CPT violation. Note that 
the measured value of sin(2/3rf) can go down from its CKM expectation for Im(<5) > 0, in 
fact, for Im(<5) « 0.07, sin(2/3 c ;) « 0.66, as can be seen from figure [51 While this value of 
Im(<5) generates a mixing phase for the B s system that is consistent with the CDF and DO 
measurements at la, one must remember that 5 need not be a flavour-blind parameter. 



5 Summary and Conclusions 

We have investigated the possibility of CPT violation in neutral B systems. CPT is a 
symmetry that is expected to be exact and the violation, even if it exists, should be quite 
small. However, it is possible to measure even a small CPT violation from the tagged and 
untagged decay rates of the neutral B mesons. In particular, for single- amplitude decay 
channels, the coefficients of the trigonometric terms sin(Amt) and cos(Amt) can effectively 
pinpoint the nature of the CPT violating parameter 5. This is an interesting possibility for 

4 We do not take the measurements coming from b — > s penguin channels because of their inherent 
uncertainties. 
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Figure 3: Variation of Aqpt with Re (5) for the system. The three lines, from top to 
bottom, are for Im(<5) = —0.1,0 and 0.1 respectively. 
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Figure 4: Variation of Acpt with Im(5) for the B,i system. The three lines, from top to 
bottom, are for Re(5) = —0.1,0 and 0.1 respectively. 
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Figure 5: Variation of sin(2/?d) with Im(5). 

the decays B s — > D+D~ and Bs —> J/'4 , 4 > (with an angular analysis). Even a small CPT 
violation, allowed by the mixing constraints for the B^ system, can make the B s mixing phase 
more compatible with the Tevatron measurements, at the level of about la. On the other 
hand CPT violation should not affect the semileptonic CP asymmetries, as the corrections 
are quadratic in nature, and expected to be negligible for small 5. Thus, a correlated study 
of the CP asymmetries in B s — > Jipcf) and B s — > D^D~ vis-a-vis B s — > D s iv might be useful 
to pinpoint the CPT violating effects. This, we feel, is something that the experimentalists 
should look for in the coming years at the LHC. 
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